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We motivate and use the concept of free random variables 
for the studies of the de-pinning transition of flux lines in su- 
perconductors as recently discussed by Hatano and Nelson. 
We derive analytical conditions for the critical points of the 
eigenvalue distribution in the complex plane for all values of 
the transverse magnetic field in one-dimension, in agreement 
with numerical calculations. We suggest a relation between 
reduced nonhermitean quantum mechanics in D-dimensions 
and nonhermitean random matrix models with weak nonher- 
miticity. 

PACS numbers : 72.15.Rn, 74.60. Ge, 05.45. +b 



1. Recently Hatano and Nelson jjj have shown that 
the de-pinning of flux lines from columnar defects in su- 
perconductors in D+l-dimcnsions, may be mapped onto 
the world lines of bosons in D-dimensions. The pinning 
and hence localization by the columnar defects is mapped 
onto an on site real randomness, and the de-pinning by 
the transverse magnetic field is mapped onto a directed 
hoping, resulting into non-hermitean quantum mechan- 
ics. While it is generally accepted that all eigenstates are 
localized in one- and two-dimensions in the presence of 
randomness, it is clear that the flux lines are de-pinned 
by a strong transverse magnetic field. The de-pinningin 
one- and two-dimensions was studied numerically in |l|. 

In this letter, we would like to show that the non- 
hermitean tight-binding model discussed by Hatano and 
Nelson in one-dimension can be analyzed in a straightfor- 
ward way using the concepts of free random variables. In 
section 2 and 3, we introduce the model and motivate the 
use of the addition law for free random variables. In sec- 
tion 4, we derive an explicit condition for the end-points 
of the distribution of (localized) eigenvalues on the real 
axis for arbitrary transverse magnetic field. In section 5, 
we extend our analysis to the complex eigenvalue plane. 
In section 6, we use the uncertainty principle to show 
that the constant mode sector of this and related mod- 
els may be amenable to non-hermitean random matrix 
models with weak nonhermiticity ||[|- 

2. Following Hatano and Nelson, we consider 
the non-hermitean tight-binding Hamiltonian in second 
quantized form for D = 1 
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where c A is a boson creation operator at site A. Through- 
out, the lattice spacing a = 1. The diagonal entries 
are random with elements distributed uniformly between 
(—A, A), and the deterministic part Hq is off-diagonal 
with hopping strengths te ±h /2 (t < 0). Here h is the 
typical strength of the 'transverse magnetic field' in units 
of the flux quantum . 
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FIG. 1. Re e - Im e for t = A = 1. 

The eigenvalues of H are complex valued for h ^ 0. 
Since H(—h) = H T (h), the left-eigenvalues are equal to 
the right-eigenvalues. Their distribution in the z-plane 
is shown in Fig. 0, for h = 0.1,0.3 and A = t = 1. 
The results are for an ensemble of 100 matrices of size 
100 x 100. For h = 0.1 the eigenvalues are mostly real 
(localized). For h = 0.3 the eigenvalues around 2 = are 
mostly complex (delocalized). For larger size matrices, 
the width of the rim shrinks to zero in agreement with 
the results discussed in [Q. 

An important question regarding the character of the 
spectrum is the occurrence of a critical hi 1 ^ for which a 
gap near z ~ sets in for the localized states. For in- 
creasing h > hi^ the eigenvalues migrate from the real 
axis to the complex plane as also discussed by Feinberg 
and Zee Q. The migration is total for h = > hi x \ 
To try to quantify this and the bulk aspects of the spec- 
trum in Fig. |l[ we will analyze ([!]) using the addition law 
for free random variables |5||3| . 

3. In this section, we proceed to motivate the use of 
free random variables for the Hamiltonian (|l]). For that, 
we note that the diagrammatic expansion for (Q) follows 
that of random matrix models to the exception that 
the V-propagators are changed to 
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with a uniform instead of Gaussian measure. Now, con- 
sider the resolvent G(z) — tiQ(z) for this model 
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where £(z) is the 1-particle irreducible (1PI) self-energy. 
Typical contributions to 'S(z) are shown in Fig. [2| They 
may be decomposed into planar connected (pi, c) and 
nonplanar 1PI ones (np). Generically, 



S(z) = ^«trVg(z)V...g(z)V) pic . 
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+ (trvg(z)v...g(z)v) np ) 

The planar piece is equal to 
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In the last equality we made use of the invariance of H 
under lattice translations. The nonplanar piece is in- 
volved. Here, we will approximate Q(z) it by its diagonal 
part G(z) ■ 1. The latter may be moved out of the trace, 
resulting into 
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This is our main approximation. 





FIG. 2. Planar and nonplanar contributions to £(z) 
dotted lines represent Q(z). 
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Having said this, we note that the calculation of the 
resolvent would be simplified considerably if we could 
show that there exists a random matrix model M. with a 
measure P(A4) such that 



(trM") pi!C = ((trV») pi , c + <trV") nj 
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Indeed, then and hence G(z) for our problem would 
simply follow from the analogue random matrix model 
Hq + M. . The main advantage is that now only planar 
graphs in the random matrix model contribute, for which 
the addition law of free random variables || apply. 

To show this, we note that all the nonplanar parts 
(trV"') can be absorbed into effective connected i-th 

\ / np 

moments. Hence, we may use the combinatorial rela- 
tions discussed in j^], between connected and ordinary 
moments for planar diagrams to reduce the equality (0) 
just to the equality of moments 



with the pertinent measures for each averaging. This is 
equivalent to finding a probability distribution P(M.) = 
cxp(— NtrV(Ai)) such that the resolvents of M. and V 
coincide. To this end we may solve the loop equation 
(technically with a single cut) of |[o) to obtain 
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A matrix model with such a potential will satisfy all our 
requirements. However to apply the methods of free ran- 
dom variables we need only to know that such a model 
exists, the explicit form of V through (JsJ) is unnecessary. 
This concludes our proof of existence. 

To summarize, our approximation scheme for (|l|) in the 
form of (|^) includes all the planar graphs and also resums 
an infinite series of (diagonal) nonplanar contributions 
for the resolvent. It is easily implemented in the matrix 
analogue using the addition law for Ho + M. . Although 
this procedure is not exact in large N, the results are close 
to the exact results for large z and large N, as shown 
in the Table for the first few moments in the Laurent 
expansion of the resolvent, that is 
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for h = 0, t = 2 and A = 1. The deviations from freeness 
start at the eight moment and are very small. 
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4. In this section we derive the resolvent G(z) along 
the real axis using the addition law for 'hermitean' matri- 
ces. 7w(A) = ImG(A + ie) is the distribution of localized 
eigenvalues. With this in mind, the resolvent for the ran- 
dom part is simply given by 
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with an inverse (Blue's function) given by By = 
AcothAz, in agreement with The resolvent for the 
deterministic part is 
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In the large N limit Gh (z) = l/y'z 2 — t 2 for z outside 
the ellipse defined by 
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and zero inside. The inverse of the resolvent (Blue's func- 
tion) is just Bh = \/l/z 2 + t 2 . 

For large z and in the planar approximation (matrix 
model), the resolvent G{z) for ([[]) along the real axis 
follows from its functional inverse B[G(z) — z through 
the 'hermitean' version of the addition law [ | , B = B Ho + 
By — 1/z. Specifically 
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and is /i-independent. Along the real axis, the end-points 
z e of the spectrum satisfy z e = B(z c ), with dB(z c )/dz = 0, 
that is 
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For A= t=l this yields z c =1.5752 or z e =1.63915. 



5. In this section, we derive analytical conditions 
for the end-points of Fig. [ij For that, we analytically 
continue G(z) to the z-plane minus the ellipse (|l3|), and 
use it to construct the potential F(z) — \ogZ/N, where 
Z(z) = (det(,2 — H)). In the planar approximation (ma- 
trix model) and to leading order in 1/N, the potential 
follows from G(z) through F = Jdz G, [Q, that is 
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The normalization Z ~ z as z — ► oo corresponds to 
G(z) -> 1/z. For z ~ 0, Z(z) = t/2e Nh (l + 0(N )). 
This behavior is related to the other branch (G(z) = 0) 
of the resolvent, for which F(z) with z ~ is constant 
and equal to h — \og2/t. Hence, a 'cusp' develops in F(z) 
where the in and out potentials cross, that is 



F[G out (z)] = ft - log 2/t = F[G in (z)} 
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The locus of the cusp (|17]) coincides with the position of 
the complex eigenvalue distribution shown in Fig. [I] as 
N — > oo. The real (localized) eigenvalues are the rem- 
nants of the hermitean addition law discussed above, and 

disappear at some critical value of h = h c ■ 

Fig. ||a shows the behavior of the critical cusp-line 
along the real axis, where the eigenvalue spectra have 
a branching point (xbranch), versus the strength h for 
which ( |l7| ) holds. At Xb ranc h(h^) — z e the localized 
states disappear from the spectrum. The dots are the 
numerically generated branching points, while the solid 
line corresponds to using Fig. |b shows the 

same along the imaginary axis for the topmost point of 
the spectrum y ma x- For large h, the resolvent drops to 
zero in the outer region like 1/z, and using ( Jl4] ) together 
with (p6h we get for y max 
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The existence of the solution for given h defines the crit- 
ical value for which the cut starts to develop in the 
inner ellipse of Fig. |l]. The dependence of he on A is 
shown in Fig. || (solid line). The dotted lines are gen- 
erated analytically by working out the leading contribu- 
tions to F in (nfl) . For small values of A (dotted line) 
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FIG. 3. The endpoints of the complex eigenvalue spectra 
versus h at A=l, t= 1. See text. 

These results are also reproduced analytically using a 
semi-circular distribution for the random part, that is 
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when r = A 2 /3 — » oo, as indicated by the thick dashed 
line in Fig. 
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FIG. 4. h c versus A numerically (solid line). For the an- 
alytical curves (dotted and dashed) see text. 

Finally, the dependence of the eigenvalue distribution 
shown in Fig. [l]on 1/N, maybe qualitatively understood 
by using the circular version of (|l]). This amounts to 
trading t/2e h ~*r and t/2e~ h — »0. For a uniform distri- 
bution of eigenvalues, the secular equation is 



Y[(X-Vi) = -(-r) N 



(22) 



For small values of A, hence Vi, (|22|) can be solved per- 
turbatively, 
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with the ansatz A, 



re 



i2nj/N 



+ e, 



giving 



-^Vi/N. Typically, (N 2 e 2 ) = N(V 2 ) = NA 2 /3, 

showing the corrections to be e ~ 1/y/N in large N. 

Numerically the width of the boundary between the 
outer (G(z) holomorphic) and the inner (G(z) = 0) 
maybe analyzed using 



SX = max (| Ai| — \Xj\) 



(24) 



along the real axis. The dependence on the size of the 
matrix is indeed ~ 1/y/N as shown in Fig. 




FIG. 5. S\ versus N (matrix size) for r — A — 1. 

6. The previous analysis borrows on some of the 
methods discussed in jjj for random matrix models. 
However, the Hatano-Nelson model differs in an impor- 
tant way from matrix models: it knows about the di- 
mensionality of space. Recently Efetov || has used su- 
persymmetric methods to argue that in 0-dimension the 
model reduces to a nonhermitean random matrix model 
with weak-nonhermiticity Q , and suggested that the re- 
duction may yield to new developments in the context of 
oriented quantum chaos. 

The reduction is actually understandable from the 
point of view of continuum quantum mechanics of con- 
stant modes. The continuum version of (|IJ) is H = 
(p + ih) 2 + V, where p is the D-dimensional momentum 
and V the random site potential. The reduction to 0- 
dimension means that V is x-indepcndent. Classically 
this would imply that p — 0. Quantum mechanically, 
however, p fluctuates. In a box of size N , p ~ W '/ y/N 
by the uncertainty principle. Hence, the reduction of H 
to iJ» in 0-dimension amounts to 
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For sufficiently random hopping, the result is a ran- 
dom matrix model with small nonhermiticity in general. 
For D=2 this is just the case of weak non-hermiticity 
discussed by Fyodorov et al. [||. For V,W chosen in 
the GOE ensemble as motivated by V real in the D- 
dimensional version of (0), our arguments suggest local- 
ization for D=2 in ( f25| ) as also noted in ||. 

The present arguments may also extend to other mod- 
els. For example, the (massless) QCD Dirac operator in a 
D-dimensional Euclidean box of volume N at finite chem- 
ical potential // is H = j D+1 (ij a \7 a + ipj D ), where 7's 
are Dirac matrices with a — 1,...,D, and V the covari- 
ant derivatives with external gauge fields. The squared 
operator, 



ff* = (iW a )(iW a ) + -<J ab [iV a , zV b ] + 2ip(iW D ) (26) 

with a ab — i[7 a ,7 b ]/2, is analogous to (p5|). For suffi- 
ciently random hopping and i\7 D ~ 1/ y/N, localization 
may take place for the GOE ensemble in D=2 dimensions 
at finite (i. Since two-dimensional QCD, in the limit of a 
large number of colors, exhibits quasi long-range order, 
this issue is worth investigating. 

7. We have shown that the Hatano-Nelson model (Q) , 
may be understood in terms of the addition law for free 
random variables. We used the latter to find an analyt- 
ical condition for the end-points of the complex eigen- 
value distributions. The results compare well with the 
numerics. The size of the spread of the complex eigen- 
value distribution is 1/y/N. We have presented generic 
arguments for how the constant modes of the model in 
D-dimensions and for sufficiently random hopping, relate 
to random matrix models with weak nonhermiticity. For 
D=2 our arguments confirm a recent observation by Efe- 
tov S. In light of this, it would be interesting to repeat 
our analysis for D = 2 for the three Wigner's ensembles. 
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